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consists  of 


straight  subsonic  edges.  The  method  is  then 
applied  to  several  examples. 


The  method  used  by  Doris  Cohen  (NACA  TN  No. 

1555)  1948)  is  applied  to  themotion  of  a  wing  in 

constant  yaw  in  the  case  that  the  planform  of  the 
only 

wing  consists  of  straight  subsonic  edges.  The 
method  is  applied  to  several  examples. 
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ROLLING  AND  YAWING  MOMENTS  FOR  SWEPT-BACK  WINGS 
IN  SIDESLIP  AT  SUPERSONIC  SPEEDS 
By  Seymour  Lampert 


SUMMARY 


By  the  use  of  methods  of  linearized  nonviscous,  conical— flow 
theory,  formulas  for  computing  the  rolling  and  yawing  moments  for  flat, 
thin,  swept^hack  wings  at  angle  of  yaw  have  heen  developed.  The  for—' 
mulas  were  derived  specifically  to  treat  that  family  of  wing  plan  forms 
for  which  all  edges  are  subsonic  (i.e.,  the  component  of  the  free-stream 
velocity  normal  to  the  edge  is  subsonic)  and  straight.  The  formulas  are 
also  applicable  to  wings  for  which  all  but  the  trailing  edge  is  sub¬ 
sonic.  The  rolling  and  yawing  moments  of  several  representative  plan 
forms  were  calculated.  All  the  wings  showed  stable  rolling  moments 
about  the  axis  of  symmetry  of  the  wing  plan  form,  and  showed  stable 
yawing  moments  about  a  vertical  axis  through  the  apex  of  the  wing. 


INTRODUCTION 


In  previous . treatments  of  the  problem  of  determining  the  rolling 
or  yawing  moment  due  to  steady  sideslip  or  yaw  at  supersonic  speeds 
(references  1,  2,  3,  and  4),  the  case  of  thin,  flat,  swept— back  wings 
with  all  edges  subsonic  (fig.  1)  —  that  is,  where  the  component  of  the 
free-atream  velocity  normal  to  any  edge  is  subsonic  —  has  not  been 
treated.  The  classes  of  swept— back  wings  studied  in  the  references 
cited  have  supersonic  trailing  edges,  where  the  Kutta  condition  is 
fulfilled  without  affecting  the  existing  flow  over  the  wing.  However, 
in  fulfilling  the  Kutta  condition  at  a  subsonic  trailing  edge  the  flow 
over  the  wing  will  be  affected.  The  Kutta  condition  may  be  realized 
in  such  cases  if  the  lift  distribution  (i.e.,  streamwlse  component  of 
the  perturbation  velocity)  is  made  to  go  continuously  to  zero  at  these 
edges.  Lagerstrom  in  reference  5  siiggested  that  this  may  be  accomplished 
by  the  superposition  of  conical-flow  fields  on  semi-infinite  flow  fields 
to  cancel  lift  distribution  beyond  the  edges.  The' methods  suggested  in 
reference  5  were  used  in  reference  6  to  calculate  the  lift  and  pitching 
moments  on  swept— back  wings  having  all  edges  subsonic.  In  a  similar 
manner  the  methods  of  reference  5  utilized  in  the  present  paper  to 
obtain  the  rolling  moments  of  such  wings  at  an  angle  of  yaw. 
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In  addition  to  satisfying  the  Kutta  condition  at  the  trailing  edge, 
it  must  also  be  satisfied  for  the  trailing  tip  in  the  case  of  the  yawed 
wing.  In  figure  2,  it  can  be  seen  that  as  the  wing  is  yawed,  one  initi¬ 
ally  streamwise  tip  becomes  a  subsonic  trailing  edge  where  the  Kutta 
condition  must  apply,  the  other  tip  becomes  a  subsonic  leading  edge  where 
the  lift  distribution  is  simileir  to  that  associated  with  the  leading  edge 
of  a  thin  flat  plate  in  subsonic  flow.  At  such  leading  edges  the  linear¬ 
ized  theory  gives  an  infinite  value  for  the  streamwise  perturbation 
velocity  (i.e.,  lift  distribution)  accompanied  by  a  horizontal  suction 
force.  These  conditions  at  the  tips  are  taken  into  accoimt  in  the  equa¬ 
tions  that  are  subsequently  derived  for  the  lift,  rolling  moment,  and 
yawing  moment. 

Although  the  formulas  in  this  paper  were  derived  for  wings  with 
initially  unraked  tips,  they  may  be  applied  to  wings  with  initially 
raked  tips  if  a  change  of  the  parameters  involving  the  slopes  of  the 
tips  is  made.  It  should  be  pointed  out  that  the  rolling  up  of  the 
vortices  around  a  streamwise  or  nearly  streamwise  edge  is  not  predicated 
by  the  linear  theory.  Therefore,  the  results  given  in  this  paper  may 
not  completely  describe  the  nature  of  the  flow  at  any  edge  which  is 
nearly  alined  with  the  free— stream  direction. 

SYMB0I5 


General 

V  free— stream  velocity 

M  free— stream  Mach  number 

3  VM2-1 

p  free— stream  density 

q  free— stream  dynamic  pressure 

^  coefficient  of  local  lift  per  unit  angle  of  attack 

Qcx 

a  angle  of  attack,  radians 

\|r  angle  of  yaw,  degrees  (positive  as  shown  in  fig.  3)  or  negative 

sideslip 
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Wing  Parameters 


A 


Co 

ct 

s 

S 

A 

X 


X6,y5 

I 


aspect  ratio 
root  chord 
tip  chord 
semispan 

area  of  ving'plan  form 

angle  of  sveep  of  leading  edge,  degrees 

taper  ratio 

Coordinates 


Cartesian  coordinates  in  the  stream  direction,  across  the  stream, 
and  in  the  vertical  direction,  respectively 

coordinates  of  the  apex  of  a  constantr-load  element  defined  by 
ray  a 

coordinates  of  the  apex  of  an  element  used  in  secondary  correc¬ 
tions 

coordinates  of  the  apex  of  the  trailing  edge 

spanwise  coordinate  referred  to  the  apex  of  a  constant-load 
element 


In  the  following,  all  slopes  are  measured  in  the  positive  sense 
counterclockwise  from  the  x  axis: 


a 

P 

X  (slope  of  ray  from  apex  of  wing)  ^ 

subscript) 

ta 

P 

X  (slope  of  ray  through  Xa,,ya)  =  P 

^^b 

P 

X  (slope  of  ray  through  x.jj,y|j)  =  p 

■^6 

P 

X  (slope  of  ray  through  Xg^y^)  =  p 

-  o  ^ 

x-Xa 

y-yb 

x-Xb 

y-ys 

X— 
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Rays  Describing  Wing  Geometry 


The  following  are  for  the  unyawed  wing  shown  in  figure  1: 

a^  P  X  (slope  of  ray  from  apex  through  intersection  of  right-span 

3s 

tip  and  trailing  edge)  - — 7 — r 

Co+(Ps/®t) 

m  P  X  (slope  of  right-span  leading  edge)  =  p  cot  A 

m^  P  X  (slope  of  right-span  trailing  edge) 

The  following  are  for  the  yawed  wing  shown  in  figure  2.  Here  the 
subscripts  1  and  2  refer  to  the  right  and  left  spans,  respectively. 


P  X  (absolute  value  of  tlie  slope  of  the  ray  from  apex  through 
the  intersection  of  the  tip  and  trailing  edge) 


mi^m2 


P  X  (absolute  value  of  the  slope  of  leading  edge) 


X  (absolute  value  of  the  slope  of  the  trailing  edge) 


In  general,  any  yawed  rays  di,  dg  which  describe  the  wing  geom¬ 
etry  may  be  obtained  from  the  unyawed  rays  d  by  the  following  rela¬ 
tions  : 


Perturbation  Velocity  Components 

■>> 

u,v,v  perturbation  velocities  in  the  x,y,z  directions,  respectively 

u^(a)  streanrvise  perturbation  velocity  for  a  triangular  ving  at  angle 

of  yaw 
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"b 


Dl 

Lo 

AL 

(AL)a 


AL* 

(AL')a 


if 


Cdl 

Cz 


value  of  t^(a)  for  wing  reflected  in  the  x  axis 

Incremental  velocity  on  any  constant-load  element 

value  of  %  for  ving  reflected  in  x  axis 

incremental  velocity  on  the.  elements  used  in  secondary  correc¬ 
tions 

value  of  u^(a)  for  a  =  —  5 

Forces  and  Moments 

drag  due  to  lift 

basic  lift  for  the  entire  wing 

correction  to  basic  lift 

correction  to  basic  lift  due  to  the  application  of  one  constant¬ 
load  element 

basic  rolling  moment  for  the  entire  wing  evaluated  about  the 
Mis  of  symmetry 

correction  to  the  basic  rolling  moment 

correction  to  the  basic  rolling  moment  due  to  the  application  of 
one  constant— load  element 

yawing  moment  about  the  apex  of  the  wing 

suction  force  per  unit  length  of  edge 

drag  coefficient 

lift  coefficient 

rolling-moment  coefficient  ^ 

yawing-moment  coefficient  ^ 

coefficient  of  yawing  moment  due  to  normal  force  (oC^) 
coefficient  of  yawing  moment  due  to  suction  force 


lift^ 

qs  J 


6 
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Elliptic  Integrals 


k  modulus  of  elliptic  integral  (value  given  in  text) 

k'  complementary  modulus  (  (value  given  in  text) 

K(k) 1  complete  elliptic  integrals  of  the  first  and  second  kinds, 
E(k) j  respectively,  of  modulus  k 

F(cp,k)\  incomplete  elliptic  integrals  of  the  first  and  second  kinds, 
E(cp,k)J  respectively,  of  moduliis  k  sind  argument  cp 

Ao(cp)  ^  [E(k)F(k',cp)  +  K(k)E(k,cp)  -  K(k)F(k’,cp)] 

ANALYSIS 


Tbe  present  report  utilizes  tlie  same  cancellation-of—lift  process 
as  was  used  in  reference  6  to  determine  the  induced  effects  due  to  the 
interaction  between  the  flows  on  the  upper  and  lower  surfaces  of  wings 
with  subsonic  tips  and  trailing  edges.  The  method  consists  of  super¬ 
posing  conical— flow  fields  on  a  semi— inf inite  flow  field  in  order  to 
satisfy  the  boundary  conditions  about  a  given  wing  plan  form.  The  analy¬ 
sis  of  the  characteristics  of  a  thin,  flat,  yawed,  swept  wing  with  sub¬ 
sonic  edges  by  this  method  falls  logically  into  two  main  divisions: 
namely,  (l)  that  of  computing  the  basic  forces  and  moments  on  a  given 
plan  form  due  to  the  flow  field  associated  with  a  semi-infinite  yawed 
lifting  triangular  wing  of  the  same  apex  angle  as  the  given  plan  form, 
and  (2)  that  of  determining  the  forces  and  moments  induced  by  the  inter¬ 
action  between  the  upper—  and  lower— surface  flows  around  the  subsonic 
tips  and  trailing  edges.  These  induced  forces  and  moments  on  the  wing 
may  be  computed  by  superposing  conical— flow  fields  beyond  the  tips  and 
trailing  edges  in  order  to  cancel  the  lift  of  the  Infinite  field. 

The  flow  fields  which  are  superposed  in  the  cancellation  process 
may  be  expressed  in  terms  of  the  perturbation  velocities  (u,v,wK  If 
Q  represents  a  perturbation  velocity  it  must  satisfy  the  linearized 
equation  for  a  nonviscous,  irrotational  compressible  fluid,  namely. 


(1-^^) 


5x^ 


=  0 


(1) 


It  can  be  shown  that  Q  may  be  expressed  as  a  function  of  a  single 
complex  variable;  the  variable  is  homogeneous  of  degree  zero,  that  is, 
constant  along  rays  from  a  fixed  point  in  space  (references  5  and  7). 

The  flow  fields  represented  by  Q  therefore  are  conical.  Several  of 
these  flow  fields  are  described  in  detail  in  reference  5.  The  boundary 
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conditions  whicla  vill  be  imposed  on  these  fields  for  purposes  of  the 
present  problem  will  be  indicated  later. 

In  order  to  simplify  the  analysis,  the  wing  plan  foi*m  has  been 
divided  into  a  number  of  separate  regions  as  shown  in  figure  2.  The 
regions  are  defined  as  follows: 

1.  The  basic  region  is  the  entire  wing  plan  form.  In  this  region 
first  approximations  to  the  lift  and  rolling  moment  about  the  axis  of 
symmetry  are  obtained  by  integrating  over  the  entire  plan  form  the  lift 
distribution  given  by  the  expression  for  a  yawed  triangular  wing  having 
the  same  apex  angle  as  the  given  plan  form.  This  lift  and  rolling 
moment  will  be  designated  as  the  basic  lift  and  basic  rolling  moment. 

2.  The  primary  regions  are  the  regions  which  are  influenced  by  the 
cancellation  of  the  basic  lift  distribution  for  a  sideslipping  triangular 
wing  beyond  the  tip  and  trailing  edges  of  the  desired  plan  form.  The 
corrections  to  the  lift  and  rolling  moment  which  occiJr  in  these  regions 
will  be  called  primary  corrections . 

3.  The  secondary  regions  are  those  which  are  influenced  by  addi¬ 
tional  corrections  which  must  be  applied  at  the  edges  of  the  wing  where 

^  -the  cancellation  flows  introduce  residual  lift  distributions  beyond  the 

confines  of  the  wing.  These  corrections  to  the  lift  and  rolling  moment 
will  be  called  secondary  corrections,  and,  in  general,  are  negligible 
■ii  for  small  angles  of  yaw. 


Basic  Lift  and  Rolling  Moment 


The  basic  lift  is  obtained  by  integrating  the  lift  distribution  of 
a  yawed  semi— infinite  triangular  wing  over  the  wing  plan  form  shown  in 
figure  2.  If  mi/p  and  ms/p  are  the  magnitudes  of  the  slopes  of  the 
rightr-  and  left-span  leading  edges  of  a  yawed  triangular  wing  with  sub¬ 
sonic  leading  edges,  the  velocity  distribution  in  the  stream  direction 
may  be  obtained  from  reference  8  as 


PoVa[  2mim2+a  ( mi-m2 )  ] 
P  v(mi— a)(m2+^ 


where 


(2) 


1 

E(ko)  [  ( l+m-L )  ( l+in.2 )  +  ] 


(3) 


8 
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and  E(ko)  is  the  complete  elliptic  integral  of  the  second  kind  of 
modulus 


k  =  2[  ( 1-mi^ )  ( l-m2^ )  ] 


The  velocity  dls  trihut  ion .  u*(a)  at  a  wing  section  normal  to  the  flow 
direction  is  shown  in  figure  3  for  typical  values  of  m  and  6.  The 
corresponding  lift  distribution  per  unit  angle  of  attack  for  the  tri¬ 
angular  wing  may  he  expressed  as 


^  ^  4u^(a) 
qa  Va 


(5) 


Integrating  equation  (5)  over  the  wing  plan  form  shown  in  figure  2 
yields  for  the  total  basic  lift  (see  appendix  A) 


qa 


2Po  s2 


4. 


I3(5)+I4(5) 

(6) 


The  values  of  Ii(5)  i=l,2,3»^  are  given  in  appendix  A  (equations 
(A8)  and  (A9)). 

The  basic  rolling  moment  about  the  axis  of  symmetry  of  the  wing  is 
given  by 


^o’  _  it  p  ,3 

qa  3  °  p2 


IIi(8) 


( 


ai^+S  / 


+  12(5)  -  113(6) 


( 


mgt-^gtV 

a2+~S  / 


14(8) 


See  appendix  A, 


equation  (A12),  for  II 1(6) 


and  113(6). 


(7) 
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.  Primary  Coirrections 


Tip  corrections  —  right  span  (region  II),—  The  following  correc¬ 
tion  applies  over  the  region  shown  in  figure  2  as  ABF,  where  AB 
corresponds  to  the  right— span  tip.  This  tip  acts  like  a  suhsohic  trail¬ 
ing  edge  since  the  streamwise  velocity  component  must  cross  this  edge  to 
leave  the  wing.  It  is  therefore  necessary  to  impose  the  Kutta  condition 
along ^his  edge.  This  condition  is  satisfied  hy  canceling  the  lift  dis¬ 
tribution  due  to  the  triangular  wing  beyond  the  confines  of  the  tip 
edge.  The  cancellation  is  accomplished  by  superposing  the  flow  fields 
of  an  Infinite  number  of  constantly  loaded  overlapping  sectors  of  infin¬ 
ite  extent  along  the  tip.  As  shown  in  figure  4,  each  of  these  sectors 
is  bounded  by  a  line  along  the  wing  tip  and  a  fixed  ray  a  =  constant 
from  the  apex  of  the  wing.  It  is  prescribed  that  each  sector  carry  a 
constant  load  proportional  to  the  incremental  change  of  the  loading  for 
the  triangular  wing  along  the  rays,  namely. 


V 

■'W 


da 


da 


(8) 


The  effect  on  the  ving  of  one  constant-load  sector,  defined  hy  a 
single  ray  a,  is  first  computed,  and  then  the  integrated  effect  of  all 
the  sectors  along  the  tip  is  determined.  The  solutions  for  constant- 
load  elements  shotild  he  such  that  they  introduce  no  change  in  the  angle 
of  attack  of  the  wing;  that  is ,  the  induced  downwash  on  the  wing  should 
he  zero.  In  the  stream  the  pressure  difference  must  vanish.  The 
boundary  conditions  for  a  given  elemental  sector  may  he  expressed  as 
follows  (fig,  4): 


1.  u  =  Ua  — 5  £  ta  ^  a 

2.  u  -  0  a  <  ta  ^  +  1 


3.  w  =  0 


-1  <  ta  <  -  S 


An  additional  condition  on  the  perturbation  velocities  is  that  they 
must  vanish  on  the  Mach  cone  from  Xg,,  ya-  If  a  general  ray  from  the 
apex  of  one  of  these  constant-load  elements  is  defined  hy 


ta  -  3 


x-^xa 


where 


10- 
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Xo 


(10) 


then  the  function  satisfying  the  conditions  specified^  given  by 
Lagers trom  in  reference  5^  niay  be  expressed  as 


u  =  r,p 


T]i(a^ta,)  - 


.  ^  cos-^  niCs-^ta) 

(a+5)(H-ta)+(ta+&)(H-a) 
(ta— a) (1—6) 


(11) 

(12) 


The  change  in  local  pressure  for  a  single  constant-load  sector  may 
therefore  be  expressed  as 


d 


- C08~^ 

rt  Va 


T|i(a,ta) 


-1  <  ta  <  -  6 


H13) 


✓ 


V 


The  Integration  of  the  lift  for  a  single  element  in  region  II  gives 


(AL)a  =  PVua 


^2+82 

3 


Gi(a) 


(14) 


where  Gi(a)  is  given  in  appendix  B  (equation  (Bll)). 

The  total  induced  lift  correction  on  this  tip  is  determined  by 
calculating  the  effect  of  canceling  with  one  sector  the  lift  density  at 
the  leading  edge  and  then  adding  back  the  Integrated  effect  of  all  the 
elements  along  the  tip  edge.  Since  the  magnitude  of  the  lift  density  at 
the  leading  edge  is  infinite^  the  quadrature  described  above  is  performed 
as  follows : 


AL  =  pVs^ 


P 


-u^(a)Gi(a)  + 


dU£;^(a) 


Gi(a)  da 


lim 

a — >  mi 


da 


(15) 
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Integrating  by  parts  and  passing  to  the  limit  gives 

w 


=  _  pvs2  Pf±^  u^(a)Gi*(a)  da  (l6) 

Similarly,  the  rolling-moment  correction  for  a  single  element  is 


(AL’)a  =  -  PV  Ua 


3 


Hi(a) 


(17) 


and  the  total  induced  moment  may  he  expressed  as 


qa 


p  Vs^(3^+S^) 
qap 


u^(a)Hi'(a)  da 


(18) 


The  derivatives  Gi’(a)  and  Hi' (a)  are  given  in  appendix  B  (equations 
(B13)  and  (B22)). 

Tip  correction  —  left  span  (region  IIl)«—  The  following  correction 
applies  over  the  region  shown  in  figure  2  as  DGC  and  includes  the  left 
tip  DC  which  may  be  considered  a  subsonic  leading  edge,  since  the 
streamwise  component  of  the  velocity  must  cross  this  edge  in  order  to 
reach  the  wing. 


In  canceling  the  lift  distribution  beyond  a  leading  tip,  it  is 
necessary  to  modify  the  previously  used  cancellation  fimction  (equation 
(ll))  to  satisfy  the  condition  of  xmiform  downwash  over  the  wing  plan 
form  (reference  5) .  The  modification  takes  the  form  of  an  additional 
term,  which  in  turn  introduces  a  singularity  in  the  lift  distribution  at 
the  leading  tip  analogous  to  that  at  a  leading  edge  in  subsonic  flow. 

In  treating  the  left-span  tip,  it  is  convenient,  as  in  appendix  A,  to 
consider  the  wing  to  be  reflected  in  the  x  axis  so  that  the  left-span 
tip  may  be  treated  in  the  first  quadrant.  In  figure  5  "the  boundary  con¬ 
ditions  for  a  single  element  and  the  regions  of  Influence  for  the  primary 
induced  corrections  are  given  for  the  reflected  wing.  The  boundary  con¬ 
ditions  are  as  follows; 


9 


1. 

II 

6  ^  ta  <  a 

2. 

O 

li 

a  <  ta  <  1 

3. 

o 

II 

> 

-  1  <  ta  <  5 

12 
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The  value  of  u^^  is  the  value  of  Uq,  for  the  wing  reflected  in  the 
X  axis.  (See  appendix  A,  equation  (Al),  for  the  required  transforma¬ 
tion.  ) 

The  ftuiction  which  satisfies  the  boundary  conditions  along  a  lead¬ 
ing  tip  is  givep  in  references  2  and  5  takes  the  form 


u  =  r.p 


uaj  r 


cos  ^ 


26 


a(l+S) 


/ 


(a-5)(l+a)(l+ta) 

^ta 


(  (19) 


where 


-  1  <  ta  <  5 


Tl2(a,ta) 


(a— 5) (l+ta)+(ta— S) (1+a) 
(ta-a)(l+5) 


(20) 


The  lift  distrihution  for  one  element  becomes 


W,  r 


V  qa y  TtVa 


cos-1  Ti2(a,ta)  - 


a(l+S) 


/ 


26  /  (a-6)(l+a)(l+ta) 


6— tj 


(21) 


The  integration  of  lift  for  a  single  sector  gives 


02+62 

(AL)a  =  pVUa,  -  G2(a) 

^  3 


(22) 


The  total  lift  correction  is  then 


AL  B^+62 

—  =  -  pV  s2  ii- — 

qa  0 


1^2 

/  ^Az( 

'^^2t 


a)G2' (a)  da 


(23) 


where  the  functions  G2(a)  and  its  derivative  G2' (a)  are  given  in 
appendix  B  (equations  (Bll)  and  (BI3),  respectively). 

The  rolling  moment  induced  hy  a  single  sector  becomes 


(AL*  )a  ~  —  S2(a) 


(24) 
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and  the  total  rolllng-ffloment  correction  takes  the  form 


•  ^ 

qo, 

The  function  H2(a)  and  its  derivative  Hg*  (a)  are  given  in 
appendix  B  (equations  (B20)  and  (B22),  respectively). 

Wake  corrections  (regions  IV  and  V).—  The  following  corrections 
apply  over  the  regions  defined  in  figure  2  as  O'EAB  and  O' Cl,  where 
O'B  and  O'C  are  both  subsonic  trailing  edges  at  which  the  Kutta  condi¬ 
tion  must  be  imposed.  Behind  these  edges,  the  lift  distribution  is  - 
reduced  to  zero  in  two  steps  (fig.  6):  (l)  canceling  lift  in  the  amount 
of  the  lift  density  at  the  apex  of  the  trailing  edge  with  a  single  semi- 
infinite  field  extending  over  the  whole  wake  region,  and  (2)  removing 
the  remaining  variation  by  the  superposition  of  infinitesimally  loaded 
elements  along  the  trailing  edges  in  the  manner  described  for  the  tip 
edges .  The  latter  elements  are  bounded  by  the  trailing  edge  and  rays 
from  the  apex  of  the  wing. 

The  correction  to  the  lift  resulting  from  step  (l)  above  for  the 
single  field  influences  both  regions  IV  nnfl  V.  If 


=  -  p 


v(p^+sg) 

qotP 


/mg 

u^j(a) 


H2’(a)  da 


(25) 


aa. 


t6  =  P  ^  (26) 

X-Xg 

refers  to  a  ray  from  the  apex  Xg,  yg  of  the  trailing  edge,  the  boundary 
conditions  which  must  be  satisfied  for  the  correction  function  are 


1. 

u  =  u^(-6) 

-  ma^  ^  tg  <  mi 

2. 

0 

II 

> 

—  mg^  >  >  mi 

These  boundary  conditions  and  the  regions  affected  by  the  corrections 
are  shown  in  figure  7*  The  function  takes  the  form 


u(t8)  =  r.p.  U5 


F(cp,k) 

K(k) 


(27) 


where 


Ik 
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cp  =  siri”^ 


and 


(l+ini^)(l+m2^) 


(28) 


(29) 


The  foregoing  expression  is  a  generalization  of  the  equation  given  in 
reference  9  for  the  symmetric  trailing— edge  correction  and  is  valid  when 
the  wing  is  yawed.  The  Ug  is  that  given  hy  the  triang\ilar-wing  solu¬ 
tion  for  a  equal  to  —  6,  or 

ug  =  u*(-6)  =  Po[2mim2-5(mi-m2)]Va 
P-v/  (mi+6)(iii2r-6) 


The  total  correction  to  the  lift  induced  hy  the  cancellation  of  the 
field  of  uniform  lift  distribution  in  the  wake  is  obtained  by  integrat¬ 
ing  equation  (27)  over  the  regions  of  the  wing  contained  within  the  Mach 
cone  from  the  apex  of  the  trailing  edge.  If  the  Mach  lines  intersect 
the  leading  edge,  as  shown  in.figiire  2,  it  is  convenient  to  integrate 
first  over  the  entire  region  O’BA*  between  the  trailing  edge  and  the 
Mach  cone  and  then  to  subtract  from  this  the  integral  over  the  region 
EAA'  between  the  leading  edge  and  the  Mach  cone.  This  procedure  is 
indicated  in  equation  (B24) .  Carrying  out  the  prescribed  operations 
gives  a  result  in  the  form 


^  ^  _  2pVugsg(p^-fS^) 
qa  qap 


I 


Roi-Rii 


s2  02+82) 


(31) 


The  values  of  Roi>  Rii/  and  R2i  are  given  in  appendix  B  (equations 
(B29),  (B32),  and  (B33)).  The  index  i=l,2  corresponds  to  the  right 
and  left  span,  respectively.  The  terms  Rqi  and  R02  constitute  the 
lift  corrections  for  the  region  boimded  by  the  Mach  cone  from  the  apex 
of  the  trailing  edge,  the  trailing  edge,  and  the  lines  which  define  the 
tip  edges.  The  additional  terms  Rii  and  R2i  are  those  required  if  the 
Mach  lines  from  the  apex  of  the  trailing  edge  cross  the  leading  edge. 

As  the  angle  of  yaw  is  increased,  the  effect  of  the  terms  Rn  £ind  R2i 
increases  while  that  of  the  terms  Ri2  and  R22  decreases. 
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The  correction  for  the  rolling  moment  about  the  axis  of  symmetry 
resulting  from  the  lift  correction  of  equation  (27)  is 


Mi-(-l)‘s  .  ,0° 

s3  (pa+BS)®  ^  s3  (p2+62)3 


(32) 


The  value  for  Rg^  is  given  in  appendix  B  (equation  (b41)).  As  in  the 
lift  correction,  the  terms  Rii,  R2I,  Rsi  apply  only  if  the  Mach  cone 
from  the  apex  of  the  ving  crosses  the  leading  edge. 

It  may  he  seen  in  figure  6  that  the  correction  effected  hy  super¬ 
posing  one  constant— load  sector  in  the  wake  cancels  the  largest  part  of 
the  triangulai>-wing  lifting  press\ire  behind  the  trailing  edge.  It  is 
necessary  to  remove  the  remaining  lift  by  superposing  additional 
constant-load  sectors  at  the  trailing  .e^e.  The  amount  of  lift  to  be 
canceled  by  these  elements  is  shown  in  figure  6  where  it  is  noted  that, 
along  the  right  span,  some  of  the  elements  will  add  to  the  lift  distri¬ 
bution  -vdiere  the  amount  canceled  by  the  single  field  is  excessive.  One 
of  the  elemental  sectors,  which  will  be  referred  to  as  "oblique"  sectors 
is  also  shown  in  figure  6. 

Considering  first  the  right-span  trailing  edge,  it  is  found  that 
the  function  necessary  to  cancel  the  lift  distribution  downstream  of 
this  edge  must  satisfy  the  following  boundary  conditions  for  each  sector 


1. 

u  =  ua 

a  < 

ta  ^ 

“it 

2. 

0 

il 

VI 

iH 

1 

ta  < 

a 

3* 

0 

II 

> 

<  ta 

<  + 

These  boundary  conditions  and  the  region  of  the  wing  affected  by  this 
correction  are  shown  in  figure  8. 

Note  that  while  the  condition  of  zero  induced  downwash  may  be 
specified  for  the  wing  area  adjoining  the  element,  it  is  not  possible  to 
satisfy  this  condition  on  the  opposite  wing  panel  if  the  Mach  cone  trace 
from  the  apex  of  the  element  should  intersect  the  opposite  trailing  edge 
However,  it  has  been  found  for  the  wings  investigated  in  this  paper  that 
the  downwash  induced  by  the  oblique  trailing— edge  corrections  was  negli¬ 
gible. 
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The  function  satisfying  the  specified  boundary  conditions  for  the 
oblique  trailing— edge  corrections  is 


UA(a) 


u  =  r.p. 


cos-1  ti„(a,ta) 


vhere 


na(a,ta)  = 


( l-a)  ( ta-fflit)-(mit-a)  ( 

(l-ffli^)(ta-a) 

^  ta  ^  1 


(33) 


m 


The  induced  lift  for  a  single  element  then  becomes 


(AL)a  =  pVUa  — - —  G3(a) 


P 


and  the  total  induced  lift  is 


qa  ^  qap  /  da 


(35) 


(36) 


vhere  Gg(a)  is  a  function  given  in  appendix  B  (equation  (Bll)). 

The  rolling-moment  correction  due  to  a  single  constant-load  sector 
is 


(AL*  )a  =  ?  VUa  s^ 


a  2  a 


H3(a) 


(37) 


and  the  total  rolling-moment  correction  becomes 


AL»  ^  V  s3(3g.>.5 

qa  ^  3  qap 


3fa2^R2\  0^1 


IJ  H  52^  ^ 

-6 


(38) 


(For  ^3(a),  see  appendix  B,  equation  (B20).) 
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’  The  corrections  for  the  lift  distribution,  rolling  moment,  and  lift 
for  the  left— span  trailing  edge  may  he  obtained  from  the  solutions  for 
region  IV  by  reflecting  the  wing  in  the  x  axis,  that  is,  replacing  8, 
and  Ua  by  —  8,  ma^,  and  Ug.^,  respectively.  The  functions  GaCa) 

and  Hg(a)  are  replaced  by  G^Ca)  and  .  The  latter  two  functions 

are  given  in  appendix  B  (equations  (Bll)  and  (B2Q)). 


Secondary  Corrections 

In  canceling  the  lift  about  a  given  8wep'l>-back-^ing  plan  form  in 
order  to  establish  the  correct  boundary  conditions,  the  primary  correc¬ 
tions  previo\isly  described  introduce  errors  that  can  influence  the  lift 
and  rolling  moment j  however,  only  in  the  case  where  the  Mach  lines  from 
the  apex  of  the  trailing  edge  cross  the  leading  edge  are  these  secondary 
corrections  of  any  considerable  magnitude.  In  this  case,  the  primary 
correction  which  may  cause  sizable  errors  is  the  initial  wake  correction 
introduced  by  the  single  field  from  the  tralling-edge  apex.  This  wake 
correction  adds  negative  lift  outboard  of  the  tip  station  aind  ahead  of 
the  leading  edge.  The  secondary  regions  affected  by  canceling  this 
negative  lift  are  shown  in  figure  2  as  EABE’  and  II’ C.  The  method  for 
making  these  secondary  corrections  to  the  lift  and  rolling  moment  is 
given  in  appendix  C. 


YAWING  MOMENT  DUE  TO  SIDESLIP 


In  order  to  compute  the  yawing  moment  for  a  thin,  flat  wing  at 
smgle  of  incidence  in  a  frictionless  flow,  it  is  necessary  to  find  the 
moment  of  the  total  force  (x— y  plane)  about  a  suitable  vertical  axis. 
For  a  swept  wing  with  subsonic  leading  edges,  the  total  force  in  the 
plane  of  the  wing  is  the  resultant  of  the  drag  force  due  to  lift  and  the 
thrust  force  due  to  the  edge  suction.  The  drag  due  to  lift  is  aL  for 
thin,  flat  wings,  €ind,  in  coefficient  form,  may  be  written  as 


CDl 


-  aCL=fl 

da 


(X2 


(39) 


where  a  is  measured  in  the  plane  of  symmetry  of  the  wing.  The  rolling 
moment  about  the  axis  of  symmetry  may  be  expressed  in  terms  of  the  yaw¬ 
ing  moment  due  to  the  drag  force  about  an  axis  normal  to  the  axis  of 
symmetry  as  follows: 


L’ 


m 


a 
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where  is  the  yawing  moment  due  to  the  drag  force  Dl.  In  this  equa*- 

tion  y  is  the  moment  arm  in  the  plane  of  the  wing  from  the  center  of 
pressure  to  the  axis  of  symmetry  of  the  wing.  In  coefficient  form^  the 
yawing  moment  due  to  the  drag  force  may  he  written  as 


C 


I 


a 


(hi) 


The  suction  forces  are  those  which  are  associated  with  the  singu- 
lairity  in  the  lift  distribution  which  occtirs  at  subsonic  leading  edges 
of  thin,  flat  wings.  In  reference  ID,  the  formula  for  finding  these 
suction  forces  is  giren  for  any  subsonic  leading  edge.  By  considering 
the  flow  two-dimensional  in  the  neighborhood  of  the  leading  edge,  it  is 
possible  to  determine  the  formula  for  the  suction  force  (per  unit  length) 
normal  to  this  edge  as 


(X 


(42) 


where  f(x)  is  the  strength  of  the  leading-edge  singiilarity  in  u^(a) 
and  m  is  p  times  the  local  slope  of  the  leading  edge. 

For  the  yawed  triangular  wing 


fi(,) .  iifaif) 
P 


1/4 


Va 


-^A“i+ni2 


E(ko)  4y/(l+mi)(l+m2)  +  (/ (1-^ni)  (l-ffla) 


^/5c  (43) 


where  i=l,2  refer  to  the  right  and  left  spans,  respectively,  and 
is  as  given  in  equation  (4).  This  value  for  fi(x)  above  is  the  same 
for  the  swept  wing  unless  the  Mach  lines  from  the  apex  of  the  trailing 
edge  intersect  the  leading  edge.  For  small  angles  of  yaw  the  effect 
introduced  by  the  Mach  lines  crossing  the  leading  edge  is  negligible  and 
not  considered  here.  The  yawing  moment  about  the  apex  due  to  the  suction 
on  one  edge  is 


N(ji  =  itp  yP  J'  ^  fi(x)^  xdx 


where  Xq^^  Is  the  chordwise  distance  to  the  wing  tip  from  the  apex. 
The  total  yawing  moment  due  to  suction  for  steady  sideslip  or-  yaw  is 
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The  integration  of  equation  (45)  yields  in  coefficient  form. 
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(45) 


•^ng 

a2 


■  ■  jL% .  r - 

6E(ko)  P  L  l+niiiiig  +  )  J  L 

^  (^j] 


The  total  yawing  moment  about  the  vertical  axis  (z  axis)  throu^  the 
apex  of  a  swept  or  triangular  wing  with  subsonic  edges  is 


Cn  _  Cnp-Cng 
a®  a2 


(47) 


It  will  be  noted  that  the  suction  effect  along  the  tip  which  becomes  a 
leading  edge  has  not  been  treated  here.  The  contribution  of  this  effect 
to  the  over— all  yawing  moment  is  small  by  comparison  with  that  of  the 
subsonic  leading  edges,  and  there  is  some  question  as  to  whether  the  tip 
suction  is  realized  in  a  real  fluid. 


APPLICATION  AND  DISCUSSION 


The  analysis  in  the  foregoing  section  has  been  directed  mainly 
toward  finding  the  rolling  moment  for  a  swept-back  wing  with  all  edges 
subsonic.  The  representative  wing  chosen  for  analysis  (shown  in  fig.  l) 
is  swept  back  63°,  has  a  taper  ratio  of  0.25,  and  an  aspect  ratio  of  4. 

In  addition,  the  rolling  moments  for  three  other  plan  forms  are  computed 
for  the  sake  of  comparison.  These  wings  have  the  same  sweepback  and 
wing  area  as  the  tapered  swept-back  wing  above  aind  are  of  relatively  low 
aspect  ratio.  These  plan  forms  are  shown  in  figure  9  and  represent  a 
triangular  or  delta  wing,  a  delta  wing  with  streamwise  tips,  and  a  so- 
called  "arrow  wing”  with  a  supersonic  trailing  edge.  In  figure  10,  the 
variation  of  rolling-moment  coefficient  with  angle  of  sideslip  is  plotted 
for  the  various  plan  forms  considered.  All  the  wings  gave  a  positive 
rolling  moment  for  a  positive  angle  of  yaw,  that  is,  the  left  tip  tended 
to  move  upward  and  the  right  tip  downward. 
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In  determining  the  rolling  moment  of  the  63^  swept  wing,  it  was 
found  necessary  to  consider  only  certain  of  the  primary  corrections  to 
the  basic  rolling  moment.  The  primary  corrections  which  contributed 
noticeably  to  the  ove3>~all  rolling  moment  include  the  tip  corrections 
and  the  wake  correction  due  to  the  single  uniformly  loaded  element.  The 
magnitudes  of  these  corrections  as  compared  with  the  total  rolling 
moment  per  degree  of  yaw  are  given  in  table  I.  The  secondary  corrections 
considered  were  found  to  be  negligible  except  for  the  correction  to  the 
primary  wake  correction  at  the  leading  and  tip  edges.  This  correction  ^ 
is  listed  in  table  I  as  secondary  effects.  In  the  case  of  the  delta 
wing  with  streaurwise  tlps^  the  only  corrections  necessary  were  the  pri¬ 
mary  tip  corrections.  For  the  arrow  and  delta  wing,  it  was  necessary 
only  to  compute  the  basic  rolling  moment.  The  rolling  moments  for  these 
two  wings  are  in  agreement  with  the  results  obtained  in  reference  1. 

In  general,  when  computing  the  rolling  moment  for  a  highly  tapered 
swept  wing  by  the  superposition  method,  it  is  necessary  to  consider  only 
the  basic  rolling  moment  and  the  primary  correction  given  by  equation 
(32).  For  the  wing  analyzed  here,  the  error  in  neglecting  all  other 
corrections  is  less  than  k  percent  for  the  range  of  yaw  angles  consid¬ 
ered.  It  should  be  pointed  out,  however,  that  the  primary  and  secondary 
tip  corrections  become  increasingly  important  as  the  taper  is  decreased. 

The  yawing-moment  coefficients  of  the  wings  analyzed  are  given  as 
a  function  of  yaw  angle  in  figure  11.  The  yawing  moments  due  to  normal 
force  and  leading-edge  suction  are  computed  about  the  vertical  axis 
through  the  apex  of  the  wing.  All  wings  exhibit  stable  yawing  moments, 
that  is,  for  a  positive  yaw  angle  the  wing  tends  to  yaw  to  the  left  or 
the  left  tip  tends  to  move  forward  and  the  right  tip  rearward.  The 
yawing  moments  for  wings  3  ^  were  in  agreement  with  the  results  of 

reference  4. 
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APPEKDIX  A 


FORMULAS  FOR  THE  INTEGRATED  BASIC  LIFT  AND  ROLLING  MOMENT 


LIFT  FOR  THE  BASIC  WING  PLAN  FORM 


The  total  basic  lift  is  obtained  by  integrating  the  lift  distribxt- 
tion  given  in  equation  (2)  over  the  plan  form  defined  as  region  I  in 
figure  2.  The  Integration  is  performed  over  four  separate  areas  shown 
in  figitre  1  as  OBA^  OBO%  OCO*,  and  OCD.  For  the  left  span  a  reflec¬ 
tion  of  the  wing  in  the  x  axis  permits  the  integration  to  be  performed 
in  the  first  quadrant.  Such  a  reflection  corresponds  to  the  following 
transformation  of  parameters: 


—6  — >  6  8-1^  3-2^  Bii  — >  mg  mg  — -$>  mi 


mi^  mg^  mg^  mi^ 


Then  the  lift  may  be  expressed  in  the  following  fonn; 


I2  =  ± 

qa  Va 


If 


aai  r.yi(a) 


aii  o 


vhere  the  Jacobian  has  the  value 


3(y,a)  ^  a‘ 


=  ^A-,(a)  =  Va) 


"Ag^^)  =  %(a)  =  %(a) 


The  limits  7i(a)  are  the  y  coordinates  of  the  intersections  of 
the  rays  from  the  apex  and  the  tip  or  trailing  edges;  yi(a)  and  y  (a 
are  at  the  trailing  edges  while  72(3)  and  74(3)  are  at  the  tip  edgei 
The  values  for  the  limits  of  integration  eire 
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yi(a)  =  Coa(mi^+6)  a)  ^^2+52] 

y2(a)  =  as  [p(a+6)]  //p^+SS 

yg(a)  =  Coa(m2^-8)  [ (ma^-a)  yp2+52]“^ 

y4(a)  =  as  [p(a-6)]  ^ 

and,  for  the  integration  with  respect  to  a. 


)(A5) 


ail  =  —  S  ai2  =  aai  -  ai^ 

ai3  =  5  ai4  =  agg  =  aa^ 

8,22  ~  ^*24  ” 

It  can  he  seen  that  the  third  and  foiirth  terms  of  equation  (A2)  may 
he  obtained  directly  from  the  first  and  second  terms  hy  applying  the 
transformation  given  in  equation  (Al) . 

The  integration  of  equation  (A2)  gives  for  the  total  basic  lift 


>  (A6) 


2^ft2  r 


■^O  „  P  +6 

—  =  2Pos2  - - 

qa  p2 


where 


(A7) 


'  /'[m2(mi  -ffli)-mi(mi  +m2)]y(mi-a)(m2+a) 

Ii(8)  =  7 - 77 - r  \  - ^ - - - 

(mi^-mi)(mi^+m2)  L 


mit(“i+“2) 


2  y(niit-“i)  (mit+®2) 


sin~i 


(mi-a)  (m2+mi^)-(m2+a)  (mi^-ffli) 

Imi^-aKmi+ma) 


a-ait 

L-. 


(A8) 


f 


90' 


♦ 


and 
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12(6)  =  — — /+  S(mi+iiig)^  (mi-a)(mg-5)H-(m2+a)(mi+5)  _ 

(mi+S)  (m2^)  L  2  y(mi+6)  (ms-b)  (a+6)(mi+m2) 


a=iiii 


Q'+S  1 

a—a^L^ 


(A9) 


The  values  for  13(6)  and  14(5)  may  be  obtained  by  applying  transfor¬ 
mation  (Al)  to  equations  (A8)  and  (A9)>  respectively. 


INTEGRATED  BASIC  ROLLING  MOMENfT 


The  basic  rolling  moment  is  obtained  by  integrating  the  moment  due 
to  the  basic  lift  distribution  about  the  axis  of  symmetry  of  the  wing 
over  region  I  of  figure  2.  The  integration  is  performed  over  the  four 
areas  of  region  I  as  in  the  case  of  the  basic  lift.  The  basic  rolling 
moment,  therefore,  is  defined  as 


V 

qa 


r-«  nSal  py^(a) 

X  /  / 

i=l  Sii  o 


^(3c>y) 

^(y^a) 


dy  da 


asi  pyAa.) 


I  r  f 

1=3 


ua7(®)  “  (a-8)  - T  dy  da 

^  a  d(y,a) 


(AlO) 


Here  the  Jacohian  has  the  same  value  as  that  given  in  equation  (A3)  and 
the  limits  of  integration  are  as  given  in  equations  (A5)  and  (a6) .  In 
the  same  way  as  in  the  lift  case,  the  third  and  fourth  terms  of  equa¬ 
tion  (AlO)  may  he  obtained  from  the  first  and  second  hy  equation  (Al) . 
The  integration  of  equation  (AlO)  yields  for  the  rolling  moment 
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Here  12(5)  and  14(6)  may  Tae  ottained  from  equation  (A9)  while  the 
value  for  IIi(6)  is 


IIi(S)  =  -  Ix(S)  - 


mi^+6 


J  /  mi _ mg  \ 

\  mi^+ms J 


mi^-a 


f  1  . 

^  1  ^  1 

1  2(m2-^iii)  1 

\mi^+m2 

(mi^-ini)(mi^+m2)  j 

«/  (mi— a)  (m2+a)  _ 


L  \mi^+m2  mi^-mi/  2  \mi^+m2  mi^-niii/ 

r  1 . 1  I  X 

L  2  \mi^+m2  mi^-ffli/JJ 


X 


*/  (mi^-^ni)  (mi^+m2) 


sin' 


-1  (mi-a)-(mi^-ini)  (m2->-a) 


mi^-a|  (mi+m2) 


a=ait 

a=— 5 
(A12) 


The  value  for  II 3  may  he  obtained  from  equation  (A12)  hy  equation  (Al) 
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APPENDIX  B 

FORMULAS  FOR  THE  PRIMARY  INDUCED  LIFT  AND 
ROLLING-^yiOMENT  CORRECTIONS 

INDUCED  LIFT  CORRECTIONS 


In  general,  th.e  equation  for  the  induced  lift  correction  at  a  tip 
or  trailing  edge  due  to  one  constant— load  sector  may  he  expressed  as 


(^)ai 


2pV 


l('ta) 


3t 


'aii 


ai 


^i 


25 

1+5 


/(a-5)(a+l)(l+ta) 
y  a2(5-ta) 


cos-1  ii^(a,ta) 


5(x,y)  - 


dy  dta 


^(x,y) 

S(y?ta) 


±  =  1,  2,  3,  ^ 

\l  i  =  2 

J 

The  Jacobian  has  the  value 


^(x,y)  ^  y-vai 
S(y>ta)  ta^ 


(B2) 


For  i=l,2  the  corrections  apply  at  the  right-  and  left-span  tips, 
respectively,  and  for  i=3,^  the  corrections  apply  at  the  right—  and 
left— span  trailing  edge.  The  limits  yi(ta)  are  the  values  of  y  at 
the  intersections  of  the  edges  of  the  wing  with  the  ray  ta/P  (figs,  k, 
5,  and  8).  For  i=2  (left— span  tip)  it  will  be  noted  that  there  are 
two  integrals  to  be  evaluated  in  equation  (Bl) .  This  added  correction 
is  required  at  any  tip  which  becomes  a  leading  edge,  as  does  the  left- 
span  leading  edge. 

For  purposes  of  integration  it  is  Convenient  to  write  equation  (Bl) 
in  the  form 
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where 


(AL)^^  =  2pV 


/■“  /■  ^ 


'Sli 


cos-1  T^^(a,ta,)  - 


2671  /(a-6)(a+l)(l+ta) 


(1+6) 


/ 


a2(&-ta) 


It, 

^(5 ,ta) 


i=l,2,3A 


5  =  y-yaj 


&i(ta)  = 


yi(ta)-yj 


/(B3) 


(b4) 


and 


ya^  =  as  [p(a+&)]”^ 

ya2  =  as  [p(a-6)]  ^  <^^+8^ 

y^  =  as  (mi  -ai^)  [^(mi^-a)  (ai^+6)  ]  ^  Vp^+B' 
ya^  =  a-s  (m2^-a2^)  [p(m2^-a)  (a2^-6)]  ^  ^p^^ 


d\i^(a)  ^^2  (a) 

=  ^2  =  Wag  =  Tia^  -  -^T- 


^  (B5) 


(b6) 


The  values  for  tji  and  1)3  can  he  obtained  from  equations  (12)  and  (3^)# 
respectively,  while  1)2  114  may  he  obtained  from  t)i  and  t)4  by 

equation  (Al) . 

The  limits  of  integration  for  equation  (B3)  are 

gi(ta)  =  s  ta  0i(a)^[p(ta-mi^)]~^  a/P+B^ 


l2(ta)  =  s  tg,  02(a)^[p(ta-®2^)]  ^  -/p^+B 
53(ta)  =  8  ta  ea(a)^[p(ta+B)]  ^ 


>(B7) 


l4(ta)  =  s  ta  04(a)  [p(ta-B)] 


— 1 


2+&2 
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are 
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gi(a)  = 


,  111  /  (a+5)(l+a) 

)  =  '  “  /  ■■ 

mi^+&  y  (iiii^+S)  (mi^+l) 


sM  = 


g4(a)  = 


- +  — 

1112^'“^  m2.^“^  LBI2 

^—a  a(iii2.j."^) 

111 

/(mi^-a)  ( l-a) 

mi^+8  a+8  a+8  ^ 

i  (1— 8)  (mi^+8) 

1  1^1 

f (ffl2t“a)(^“a) 

i  t  ^  ^  \  e 

(ar-S)(l-fa) 

12^-6)  (ma^+l) 


ar-8 


(1—8) (m2^— 8) 


>  (B12) 


For  the  total  induced  lift  (equations  (l6)  and  (23))  the  derivatives  of 
Gj[(a)  i=l,2  are  also  required.  These  are  obtained  hy  setting 


Gi*(a)  =  2ei(a)ei»(a)gi(a)  +gi*(a)ei(a)^ 
where  the  derivatives  0^’ (a)  are 

mi+.+8 

ni2j.— 6 

and  the  derivatives  gj^’  (a)  are 

1  r  (l+a)(mxt+5)+(a+8)(mit+l) 


(B13) 


y  (BlU) 


g2’(a)  = 


(mi^-a)2  1 

L  2  */  (init+8)  ( 1+mi^)  (a+8)  ( l+a) 

1 

-  (l+a)  (iii2^— 8)+(a— 8)  (m2^+l) 

(me^-a)^ 

^2k/ (m2.j.— 8)  (l4-m2^)  (a— 6)  (l+a) 

8 

.  (a-8)-8(a+l). -  . 

-  1 


;  (B15) 


2a^(ni2^~8)  */  (iii2^~8)  (l+ms.^)  (a— 8)  (l+a) 
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nroUCED  ROLLING-WOMENT  CORRECTIONS 


The.  induced  rolling-ffloment  correction  for  a  single  elemental  sector 
is  defined  as  the  moment  of  the  induced  lift  of  one  element  about  the 
axis  of  symmetry  of  the  wing.  The  general  equation  for  the  rolling 
moment  for  one  elemental  sector  is,  therefore. 


jr 


'aii  7a± 


cos-1  T]^(a,ta)  - 


2Syi 

(1+S) 


.  /(ar-5)(lfa)(n.tl)  1  /  ^ 


ta) 


i=l,2,3A 


As  in  the  case  of  the  induced  lift,  the  terms  for  i=l,2  are  the 
induced  moments  for  elemental  sectors  at  the  right-  and  left-wing  tips, 
while  the  terms  for  i=3^^  are  the  induced  moments  for  the  right  nnri 
left  trailing  edges.  The  values  of  ^(x,y)/d(y,ta)  and  ua^,  are  given 
in  equations  (B2)  and  (b6),  respectively. 


For  purposes  of  integration,  equation  (B16)  may  be  put  in  more 
convenient  form;  namely. 


(AL*  )ai 


/  (a-^)  (a+1)  (1-i-ta) 

Ti+^  J  a2(B-ta) 


Tii(a,ta)  — 

(■ta)+yaif2(ta)] 


^  d|  dtaXBlT) 


i=l,2,3A  j 


where  |  and 
are  defined  as 


are  given  in  equation  (b4)  while  fi(ta)  and  fsCta) 


t„+5 


f2(ta)  = 


a+S 

a 


(B18) 
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The  limits  of  integration  are  given  in  equations  (B7)  and  (b8)  while  the 
values  for  ya^  are  given  in  equation  (B5).  The  integration  of  equa¬ 
tion  (BI7)  yields  for  the  induced  rolling  moment  of  one  element  at  a  tip 
or  trailing  edge 


(AL'  )  aj^ 


pVai  s3 


Hi(a) 


i=l,2,3,i^ 


Here  the  functions  Hi(a)  are  defined  as 

Hi(a)  =  ei(a)\i(a) 
i=l,2,3,4 


(B20) 


where  the  functions  0j(a)  are  defined  in  eqxaation  (B9)>  while  the 
functions  hi(a)  are  given  as  follows; 


hi(a)  =  gi(a) 


a+&  1  1-6 


1  +  Qi(a) 

2  V  a  2  1+mi^ 


)]- 


Ql(8’)  (^~S)  (^+5) _ 

6(l+mi^) (mi^+6) (mi^-a) 

OgCa)  (l+S)  (a-^) 


+ ^  g2i(»)  7^2^ 

3  (mat-a) 


6(l+m2^) (m2^— 6) (m2^— a) 


> (B21) 


g2i(a)  - 


113(a)  = 


h4(a)  = 


III2 

(mii— ait) (a+8) 

2  +  ■  ■  . . . -  u 

(mi^— a) (ai^+8) 

(®2t~®’2t)  (s-~S) 

2  +  . . ■■■■—  ■■  ■  I 

(m2^— a)  (a2^— 6) 


(a— 6) (l+a) 
(m2^-8)(l+m2^) 

g3(a) 

g4(a) 
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The  values  for  g^Ca)  and  (a)  are  given  in  equations  (B12)  and  (B15), 
respectively. 


IMDUCED  LIFT  CORRECTION  IN  WAKE  DUE  TO  THE 
SINGLE  CONSTANT-LOAD  ELEMENT 


The  total  lift  correction  In  the  wake  due  to  the  single  constant¬ 
load  element  is  given  as 

^  =  V  r  /  45  it, _ 


AL  _  2pVUf.  ^  P  f' 


^i,%) 


i=l 


.^i2 


f  f  rC^.k)  «  its 

Jtst  All 


where 


(B24) 


X-Xg  X-Xg 


(B25) 


and  the  Jacohian  has  the  value 


^(x,S)  ^  _Pi_ 


(B26) 


The  limits  ts^^  refer  to  the  rays  from  the  apex  of  the  trailing  edge 
through  the  leading— edge  tip.  For  the  right  span  t5j^  has  the  value 


_  smi  ( ) +CoSp (mi+5 ) 

'  s(p2+52)-CoP(mi+&) 

while  the  limits  for  the  integration  with  respect  to  £  are 

^  ^  Cots(mi+5) 

( )  </ p^+5® 

^  stsV^g+52 
3(t6+8) 


(B27a) 


(B27b) 
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The  limits  5 21^  a^<3.  lag  for  the  left  span  may  he  obtained  hy 

applying  eq\iation  (Al)  to  eqiiation  (B27). 

O’ 

Performing  the  integration  indicated  in  equation  (B24)  gives 


qa 


2PVU6  ^  saQW) 
qa  p 


I 


Where  the  values  of  Roi^  Eii^  'and  R2I  are  as  follows; 


where 


Roi 


(1— 8)(mi^+6) 


Jt  Ao((po)  j  (Hmi^)(l4-iiig^) 

2  K(k)  J  (l— 5^)  (mi^+&)  (m2^— 8) 


(B29) 


^o^'Po)  =  I  t^(«PoA’  )E(k)  -  F(q)o,k»  )K(k)  +  E(cpo,k’  )K(k)]  (B30) 


The  values  of  Aq  for  various  values  of  cp^  are  tabulated  in  refer¬ 
ence  11.  The  value  for  R02  may  be  obtained  from  equation  (B29)  by 
making  the  transformation  of  parameters  given  in  equation  (Al) 


F((PiA)  r  l+tSi  1  /  (l+mi^)(l+m2^) 

K(k)  ^  (l— 8)  (tg^+S)  K(k)  ^  (l~8^)  (mi^+&)  (m2^— 6) 


(B32) 
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R21  - 


vhere 


^  K(k)  ^ 


l+tSi 


(1+mi) (tB^-mi) 


(l+mit) (l+®2t) 


K(k)  y  (I^iii2)(mi^-^ni)(iii2^+mi) 


F(cPiA)J 


(B33) 


=  tan" 


/  2  !•  (jfi  sin  ii|a<  sIeIl  nv^ 
-1  ^=.1 _ ! _ 1. 


1=1,2 


2 

1+2  E  cj^  cos  njJi^  cosh. 
n=l 


y  (B3l^) 


—  «■ 


K(k') 


CD  =  e 


(B35) 


til  =  ^2  =  « 


V,  =  rt 


Vo  =  Jt 


K(k) 

F(cPo,k») 

K(k) 

F(q)2ik') 

K(k) 


(B36) 


}  (B37) 


cpi  =  sin 


— 1 


(H-mit)(l-t.&i) 

( 1-nii )  (nii.t;+m2t ) 


(B38) 


<P2  =  sin 


— 1 


(m2f*-mi)(l+mit) 

(l+mi)(mi.j.+m2.t.) 


(B39) 
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The  values  for  o)  are  tabulated  in  reference  12,  while  k,  k’,  cpQ,  and 
tgj^  are  given  in  equations  (B3I)  and  (B27).  It  should  he  noted  that 
when  cPi  equals  rt/2,  Ji  is  zero.  The  values  for  R12  and  R22  may 
he  obtained  from  equations  (B32)  and  (B33)>  respectively,  if  equation 
(Al)  is  applied,  .  . 


IKDUCED  R0LLIRG-410MENT  CORRECTION  AT  ROOT  CHORD 


The  induced  rolling-moment  correction  corresponding  to  the  lift 
correction  in  eqtiation  (B28)  is 


AL* 

qa 


2  Uss3(pg+6g) 

3  ^  qap 


X  { 


Roi  -  Rii  + 


i=l 


Co®  [mi-(-l)i6]  „3  „  .  Co®  [mi-(-l)i6] 


s®  (p2+52) 


P  R2i  + 


s3  (p2+62)' 


p®  R3i|(BkO) 


The  expression  Rai  is 


^  (mi+5)^  R(q)i>^)  J  r  (1+mi)^  -  1  + 

(l-^Hi^)  K(k)  ^  ^  (hSj— mi)®  (l-mi)  (mi^-HUi) 


(1-mi) 


^  -  J. 


J'(<I>l,k)  J 


3  -  2  mi  + 


(I-T2)  (T2-k2)  J 


/ 


(l-T2)(T2Ht  )  (l-T2)(T2-k  )E(Cpi,k)  L 


E((pi,k)  - 


TgTi  V^(l-Ti^)(l-k^T?y 

1_T2Ti2 


} 


(b41) 


where 
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Tl  = 


T2  = 


J  ( ) 


y (b42) 


The  form  for  Bsa  may  he  obtained  by  using  trsins  format  ion  (Al) . 
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APPENDIX  C 

FORMULAS  FOR  THE  SECONDARY  INDUCED  LIFT  AND 
ROLLING-MOMENT  CORRECTIONS 


The  corrections  discussed  here  are  those  vhich  are  made  at  the 
leading  and  tip  edges  to  correct  for  the  residual  lift  and  rolling 
moment  introduced  hy  the  single  constant— load  correction  in  the  wake. 

The  formula  for  the  streamwise  velocity  for  a  given  elemental 
sector  at  the  leading  edge  of  the  wing  is  of  the  same  form  as  eqxia— 
tion  (19),  and  for  the  right  span  is 


QQg-i  (ts-^ni)  (l+tb)+(th-mi)  (1+tfi)  ^  2mi</(tsr-ffli)(l+t6)  /l+tb  1 

(tbr-t6)(l-Hai)  t5(l+mi)  J  mi— tfa  J 

(Cl) 


where,  if  y^  is  the  apex  of  the  element. 


tb  =  P 


y-yp 

X-Xb 


(C2) 


and 


tg  = 


*h-*6 


(C3) 


^h 


d 

dtg 


F((p,k) 
K(k)  . 


dtg 


(C4) 


The  term  in  brackets  to  he  differentiated  is  uCtg)  of  equation  (2l^)i 
The  equivalent  left— span  solution  may  he  obtained  from  equation  (Cl)  hy 
applying  the  transformation  of  equation  (Al) .  In  order  to  apply  equa¬ 
tion  (Cl)  at  the  right  tip,  it  is  necessary  to  replace  mi  hy  — 6  «nH 
at  the  left  tip  hy  8.  It  will  he  noted  that  at  the  right  tip  the 
second  term  of  equation  (Cl)  is  imaginary,  while  at  the  left  tip  it  is 
real.  Therefore,  it  is  necessary  to  retain  the  second  term  in  the  square 
brackets  for  evaluating  the  correction  along  the  left  span. 
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Induced-Lift  Corrections 


The  secondary  correction  to  the  lift  for  a  single  element  at  the 
leading  edge  of  the  right  span  takes  the  form 


dti 


(AL), 


L.E. 


[cos  ■ 


Tl(ti3,mi)  - 


% 

u^  [cos~^  Ti(t-b,mi)  - 

Ti(ti3,mi)]  d|  dti3  I 


(C5) 


where  In  and  I21  correspond  to  the  intersection  of  the  rays  t-jj  with 
the  trailing  and  tip  edges,  respectively,  of  the  wing  plan  form,  and 
tb^  is  the  ray  tb  passing  through  the  trailing-edge  tip.  The  expres¬ 
sion  in  "brackets  denotes  the  "bracketed  expression  in  equation  (Cl) .  The 
limits  of  integration  have  the  values 


6  _  Co'^)(°‘i'*'8)  ("ts-iPi-t.) 

^  ^  i/p2+52  (mi^-t-b )  ( t6-mi ) 

tj,  [s(p^+6^)(t5-ffli)-fCo(t5+S)(mi->-5)] 

3  (tb+6)(t5-ffli)  7^2+52 

P  Co"t5  ( nii+6 )  ( )  ( "tg^i ) 
pco(mi+6)  (mi^+6)-s(32+6^)  (t6-^i) 


(C6) 

(C7) 

(C8) 


while  the  left-span  correction  may  be  obtained  from  equation  (C5)  Ly 
applying  equation  (Al) . 

The  secondary  correction  to  the  lift  along  the  right— span  tip  for  a 
single  element  is 


MCA  Tn.2262 


39 


(AL)i,.  dtg  = 


Uv,  [cos  ^  il(t-b,-8)  - 


T(tt,-&)]  d|  dt^, 


The  limit  of  Integration  l^i  refers  to  the  intersection  of  the  ray  t-jj 
yith  the  trailing  edge,  and  has  the  value 


st-b(tg-mi^)  </pg+Sg 
p(t5+6)(mit-tT3) 


(CIO) 


The  leftr-span-tip  correction  may  he  obtained  hy  applying  equation  (Al) 
to  eqiiation  (C9) . 

The  total  secondary  correction  to  the  lift  is  obtained  by  summing 
the  lift  due  to  all  the  elemental  sectors  along  the  leading  and  tip 
edges.  The  correction  along  the  tip  includes  an  initial  sector  of  con-^ 
stant  magnitude  F[<p(t5j^),k/K(k)]u5.  It  should  be  pointed  out  that  due 
to  the  different  orientation  of  the  elemental  sectors  along  the  tip  and 
leading  edges  there  is  an  overlapping  of  sectors  in  the  region  bounded 
by  the  extensions  of  the  leading  edge  and  the  ray  t5 from  the  apex  of 
the  trailing  edge  through  the  tip  of  the  leading  edge.  An  approximate 
method  of  correction  for  this  overlapping  is  to  let  the  Initial  tip 
sector  be  bounded  by  the  extension  of  the  leading  edge  rather  than  by 
tSj^.  Then  the  total  secondary  correction  to  the  lift  for  the  right  span 
is  obtained  as 


3-  (al)  dtg  +  '^*8  + 


2pVP  ^  F[T(t6j^)>k]  p  ®  r^^^fcos-i  (mi+5) (l+tT3)+(ti3+S) (l+mi) 
Jt  K(k)  Jo  L 


mi(l-8) 


-(S+tb) 


^  d|  dtb 


(Cll) 


where 


.  ^  mis(p^+5^)  +  Cr,5p(mi+5) 

s(p^+8^)—  CoP(iiii+8) 


(C12) 
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The  total  secondary  correction  to  the  lift  at  the  left  span  may  he 
obtained  from  equation  (Cll)  by  transformation  (A5) .  If  the  Mach  cone 
from  the  apex  does  not  cross  the  leading  edge,  then  it  is  necessary  to 
evaluate  only  the  second  Integral  of  equation  (Cll),  since  the  limit 
tg^  becomes  unity. 


Induced  Rolling  Moment 


The  secondary  correction  to  the  rolling  moment  for  an  elemental 
sector  at  the  leading  edge  is 


dtg  ^*5  “  ;7p;6'g{/  ‘V  - 


1  0 


Ti(tb,mi)]  t6f;^^(tb)+(S6^)fg^(tb)]  g  ^ 

tb^ 


pTUx  nisi 

J  J  f  cos-1  ,^(tb,mi)  -  Tj^(tb,mi)]  X 

■fcbi  o 

Ufii(tb)+(S  5,  )fgi(t'b)] 


^^b‘ 


S  dg  dtb 


} 


(CI3) 


where  fii(tb)  and  f2i(tb)  are  defined  as 


fii  = 


fai  - 


t-b+s 

tb 


*6 


} (ci4) 


and 


I  =  °o  ts(mi+S) 

tJW+^  tg-mi 


(CI5) 


The  limits  of  integration  are  as  given  for  the  lift  in  equations  (C6), 
(C7),  and  (c8) . 
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Th.e  secondary  correction  to  the  rolling  moment  for  a  single  element 
along  the  right  tip  takes  the  form 


2pVB' 


“(AL'),  at8  =  , -  , 

^tip 


/— o  pSai 

/  xijj  [cos-1  Ti(th,-8)- T(t|j,-5)]x 


[  5f ii(th)+( 56i)fsi(th) ] 

V 


I  d|  dtb 


The  value  for 


in  equation  (Cl6)  is 


stg  V  p2+82 
=  p(t8+8) 


(C16) 


(CIT) 


Tile  left— span  correction  is  obtained  from  equation  (Cl6)  by  transforma¬ 
tion  equation  (Al) .  The  limit  of  integration  is  the  same  as  pre¬ 
viously  given  In  equation  (CIO). 

The  total  secondary  rolling-moment  correction  for  the  right  span  is 


AL^ 

qa 


-I  ^  (AL»)i  dt5  +  r  ^  (AL’)i  dts  + 
ia\  /  dtg  ^L.E.  J  ^  ots  tip 

1  *5i 

F[cp(t8i),k] 


t 


qa  ^ 

2pVP^ 

n 


K(k) 


/  / 


C(mi,5,tb)  X 


[6fii(th)  +  (I8i)f2i(th)3  <il  ith  I 


(CI8) 


where  the  limit  tg^  is  as  given  for  the  lift  in  eqxiation  (C12),  the 
value  for  is  given  in  eqtiation  (C17)>  aud  c(mi,5,tt))  is  the 

bracketed  expression  in  the  last  integral  in  equation  (Cll) .  The  left- 
span  correction  is  again  found  by  applying  equation  (Al) .  If  the  Mach 
cone  from  the  apex  of  the  trailing  edge  does  not  cross  the  leading  edge 
then  all  terms  of  equation  (CI8)  but  the  second  vanish. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Oct.  30,  I95O. 
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Basic  Region 
Region  I  OABO'CD 
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Figure  2~  Regions  over  which  the  various  corrections  apply. 
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Figure  3.-  Velocity  distribution  u^(a)  for  a  yawed  triangular  win 
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of  attack  with  yaw  angle  v  for  M-V? 
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Yaw  angle  in  degees 


Figure  n.— Variation  in  yawing-moment  coefficient  per  unit  angle 
of  attack  with  yaw  angle  for 
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